Non-linear QCD dynamics and exclusive production in ep collisions 
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The exclusive processes in electron-proton (ep) interactions are an important tool to investigate 
the QCD dynamics at high energies as they are in general driven by the gluon content of proton which 
is strongly subject to parton saturation efltects. In this paper we compute the cross sections for the 
exclusive vector meson production as well as the deeply virtual Compton scattering (DVCS) relying 
on the color dipole approach and considering the numerical solution of the Balitsky-Kovchegov 
equation including running coupling corrections. We show that the small- 2: evolution given by this 
evolution equation is able to describe the DESY-HERA data and is relevant for the physics of the 
exclusive observables in future electron-proton colliders and in photoproduction processes to be 
measured in coherent interactions at the LHC. 
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I. INTRODUCTION 

The understanding of the high energy (small x) regime 
of Quantum Chromodynamics (QCD) has been one of the 
main challenges of this theory, which has been intensely 
investigated through high energy collision experiments. 
This regime, where one expects to observe the non-linear 
behavior predicted by theoretical developments, has been 
explored in ep collisions at DESY-HERA and pp/dA col- 
lisions at BNL-RHIC and, in a near future, mpp/pA/AA 
collisions at CERN-LHC. In particular, at high energies, 
the growth of the parton distribution is expected to sat- 
urate, forming a Color Glass Condensate (CGC), whose 
evolution with energy is described by an infinite hier- 
archy of coupled equations for the correlators of Wil- 
son lines [H [1]. In the mean field approximation, the 
first equation of this hierarchy decouples and boils down 
to a single non-linear integro-differential equation: the 
Balitsky-Kovchegov (BK) equation [l|, This equa- 
tion determines, in the large- iVc (the number of colors) 
limit, the evolution of the two-point correlation function, 
which corresponds to the scattering amplitude M{x, r, h) 
of a dipole off the CGC, where r is the dipole size and 
h the impact parameter. This quantity encodes the in- 
formation about the hadronic scattering and then about 
the non-linear and quantum effects in the hadron wave 
function (For recent reviews, see e.g. Q). Recently, the 
next-to-leading order corrections to BK equation were 
calculated [^-Q through the ressumation of UgNf contri- 
butions to all orders, where Nf is the number of flavors. 
Such calculation allows one to estimate the soft gluon 
emission and running coupling corrections to the evolu- 
tion kernel and, in particular, the authors have verified 
that the dominant contributions come from the running 
coupling corrections, which allows to determine the scale 
of the running coupling in the kernel. The solution of 
the improved BK equation was studied in detail in Refs. 
[1, Q . Basically, one has that the running of the coupling 
reduces the speed of the evolution to values compati- 



ble with experimental data, with the geometric scaling 
regime being reached only at ultra-high energies. In Ref. 
, the solution of the improved BK equation was used to 
calculate the pseudorapidity density of charged particles 
produced in nucleus-nucleus collisions and a remarkable 
good agreement with the RHIC data was observed. More 
recently, a global analysis of the small x data for the pro- 
ton structure function using the irnproved BK equation 
was performed [1] (See also Ref. [l^). In contrast to 
the BK equation at leading logarithmic asln(l/a;) ap- 
proximation, which fails to describe data, the inclusion 
of running coupling effects to evolution renders BK equa- 
tion compatible with them. The improved BK equation 
has been shown to be really successful when applied to 
the description of the ep HERA data for the inclusive 
and diffractive proton structure function d-IHl, as well 
as for the forward hadron spectra in pp and dA colli- 
sions [ll|, , which motivates us to extend the study for 
exclusive observables. 

Exclusive processes in deep inelastic scattering (DIS) 
have appeared as key reactions to trigger the generic 
mechanism of diffractive scattering. In particular, the 
diffractive vector meson production and deeply virtual 
Compton scattering (DVCS) have been extensively stud- 
ied at HERA and provide a valuable probe of the QCD 
dynamics at high energies. In a general way, these pro- 
cesses are driven by the gluon content of target (proton 
or nuclei) which is strongly subject to parton saturation 
effects as well as considerable nuclear shadowing correc- 
tions when one considers scattering on nuclei. In par- 
ticular, the cross section for exclusive processes in DIS 
are proportional to the square of scattering amplitude, 
which turn it strongly sensitive to the underlying QCD 
dynamics. They have been successfully described us- 
ing color dipole approach and phenomenological model 
inspired in general aspects of parton saturation physics 
13- 15]. Here, we will make use of numerical solution of 
the Balitsky-Kovchegov equation including running cou- 
pling corrections in order to estimate the contribution of 
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the saturation physics for exclusive processes. Our anal- 
ysis is relevant for the physics to be studied in future 
electron - proton collider, as e.g. the LHeC and 
in photoproduction processes in coherent interactions at 
the LHC [ITJ. This paper is organized as follows. In next 
section (Section |lTl) the main formula for computing the 
differential cross section for exclusive processes in DIS 
are presented. Moreover, we discuss the main aspects of 
the running corrections for the BK equation. In Section 
mil we present our results and discussions. 



II. EXCLUSIVE PROCESSES IN DIS AND THE 
RC BK SOLUTION 



Let us consider photon-hadron scattering in the dipole 
frame, in which most of the energy is carried by the 
hadron, while the photon has just enough energy to disso- 
ciate into a quark-antiquark pair before the scattering. In 
this representation the probing projectile fluctuates into 
a quark-antiquark pair (a dipole) with transverse separa- 
tion r long after the interaction, which then scatters off 
the hadron . In the dipole picture the amplitude for 
production of an exclusive final state E, such as a vector 
meson {E — V) oi a. real photon in DVCS {E = 7) is 
given by (See e.g. Refs. 

AfP^^'P (x, g2, A)= Jdz d'r i^^*^)T,L A,^ix, r, A) ,(i; 

where (^'■^*^')t.l denotes the overlap of the photon and 
exclusive final state wave functions. The variable z (1 — z) 
is the longitudinal momentum fractions of the quark (an- 
tiquark) , A denotes the transverse momentum lost by the 
outgoing proton (t — — A^) and x is the Bjorken variable. 
For DVCS, the amplitude involves a sum over quark fla- 
vors. Moreover, Aqq is the elementary elastic amplitude 
for the scattering of a dipole of size r on the target. It 
is directly related to J\f{x, r, b) and consequently to the 
QCD dynamics (see below). One has that [ij 



where /3 is the ratio of real to imaginary parts of the scat- 
tering amplitude. For the case of heavy mesons, skewness 
corrections are quite important and theyare also taken 
into account. (For details, see Refs. [l^.Tl9| and Section 

The photon wavefunctions appearing in Eq. ([3]) are 
well known in literature |14|. For the meson wavefunc- 
tion, we have considered the Gauss-LC model [14] which 
is a simplification of the DGKP wavefunctions. The mo- 
tivation for this choice is its simplicity and the fact that 
the results are not sensitive to a different model. In pho- 
toproduction, this leads only to an uncertainty of a few 
percents in overall normalization. We consider the quark 
masses mu^d,s = 0.14 GeV, rric — 1.4 GeV and mt, — 4.5 
GeV. The parameters for the meson wavefunction can 
be found in Ref. ^^]. In the DVCS case, as one has a 
real photon at the initial state, only the transversely po- 
larized overlap function contributes to the cross section. 
Summed over the quark helicities, for a given quark flavor 
/ it is given by 21], 



Aqqix,r,A) = i J d'^b 



^^e-'^-^2N{x,r,b) , (2) 



where b is the transverse distance from the center of the 
target to one of the qq pair of the dipole. Consequently, 
one can express the amplitude for the exclusive produc- 
tion of a final state E as follows 



A 



T.L yx.Q = i 

X (*^^')t 2A/'(a;, r, 5 



(3) 



where the factor — z)r].A in the exponential arises 
when one takes into account non-forward corrections to 
the wave functions [20|. Finally, the differential cross 
section for exclusive production is given by 

^ij*p^Ep) = ^|<r''^(.,Q^A)p(l-|-/32),(4) 



Nrap. 



27r2 

mjKoieir)Ko{e2r)} , 



(5) 



where we have defined the quantities 612 — ^^ Qi 2 + 
and z = (1 — z). Accordingly, the photon virtualities 
are Qf ~ (incoming virtual photon) and Q2 — 
(outgoing real photon). 

The scattering amplitude J\f{x,r,b) contains all in- 
formation about the target and the strong interaction 

S' lysics. In the Color Glass Condensate (CGC) formalism 
y] , it encodes all the information about the non-linear 
and quantum effects in the hadron wave function. It can 
be obtained by solving an appropriate evolution equation 
in the rapidity y = ln(l/a;), which in its simplest form is 
the Balitsky-Kovchegov equation. In leading order (LO), 
and in the translational invariance approximation — in 
which the scattering amplitude does not depend on the 
collision impact parameter — it reads 



dJVjr, Y) 
dY 



= / dr 1 if (r, r 1 , r2 ) [AA(ri , r ) + AA(r2 , r ) 



Af{r,Y) ~ Af{n,Y)Af{r2,Y)], 



(6) 



where A/'(r, Y) is the scattering amplitude for a dipole (a 
quark-antiquark pair) off a target, with transverse size 
r = |r|, Y = lii(xo/x) (xq is the value of x where the 
evolution starts), and r2 = r — ri. K^'~' is the evolution 
kernel, given by 



^^''(r,ri,r2) 



27r^ 



(7) 
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where ag is the (fixed) strong coupling constant. This 
equation is a generalization of the linear BFKL equa- 
tion (which corresponds of the first three terms), with 
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FIG. 1: (Color online) Energy dependence of the 7p cross sec- 
tion for production for different photon virtualities. Data 
from (a) ZEUS and (b) HI collaborations lii,!!^!. 



FIG. 2: (Color online) Energy dependence of the 7p cross 
section for J/^' production for different photon virtualities. 
Data from (a) ZEUS and (b) HI collaborations IHjH^. 



the inclusion of the (non-linear) quadratic term, which 
damps the indefinite growth of the amplitude with en- 
ergy predicted by BFKL evolution. It has been shown 
[221 to be in the same universality class of the Fisher- 
Kolmogorov-Pertovsky-Piscounov (FKPP) equation ^] 
and, as a consequence, it admits the so-called traveling 
wave solutions. This means that, at asymptotic rapidi- 
ties, the scattering amplitude is a wavefront which trav- 
els to larger values of r as y increases, keeping its shape 
unchanged. Thus, in such asymptotic regime, instead 
of depending separately on r and Y, the amplitude de- 
pends on the combined variable rQs{Y), where QsiY) is 
the saturation scale. This property of the solution of BK 
equation is a natural explanation to the geometric scal- 
ing, a phenomenological feature observed at the DESY 
ep collider HERA, in the measurements of inclusive and 
exclusive processes f23 - [27j . Although having its prop- 
erties been intensely studied and understood, both nu- 
merically and analytically, the LO BK equation presents 
some difficulties when applied to study DIS small-x data. 
In particular, some studies concerning this equation [23- 
l32j have shown that the resulting saturation scale grows 
much faster with increasing energy {Ql x~^, with 
A ~ 4.88Ncas/TT ~ 0.5 for as = 0.2) than that extracted 
from phenomenology (A ^ 0.2—0.3). This difficulty could 
be solved by considering smaller values of the strong cou- 
pling constant as , but this procedure would lead to physi- 
cally unrealistic values. One can conclude that higher or- 
der corrections to LO BK equation should be taken into 
account to make it able to describe the available small-x 
data. 

The calculation of the running coupling corrections to 
BK evolution kernel was explicitly performed in 01 , 
where the authors included UgNf corrections to the ker- 
nel to all orders. The improved BK equation is given in 
terms of a running coupling and a subtraction term, with 
the latter accounting for conformal, non running coupling 
contributions. In the prescription proposed by Balitsky 



in |7| to single out the ultra-violet divergent contributions 
from the finite ones that originate after the resummation 
of quark loops, the contribution of the subtraction term 
is minimized at large energies. In 9] this contribution 
was disregarded, and the improved BK equation was nu- 
merically solved replacing the leading order kernel in Eq. 
^ by the modified kernel which includes the running 
coupling corrections and is given by Q 



i^^"'(r,ri,r2) 



27r2 



1 f c^sjrl) 



1 



(8) 



From a recent numerical study of the improved BK equa- 
tion 0, it has been confirmed that the running coupling 
corrections lead to a considerable increase in the anoma- 
lous dimension and to a slow-down of the evolution speed, 
which implies, for example, a slower growth of the satura- 
tion scale with energy, in contrast with the faster growth 
predicted by the LO BK equation. Moreover, as shown 
in [9, 11, 12] the improved BK equation has been shown 
to be really successful when applied to the description of 
the ep HERA data for the inclusive and diffractive pro- 
ton structure function, as well as for the forward hadron 
spectra in pp and dA collisions. It is important to em- 
phasize that the impact parameter dependence was not 
taken into account in Ref. @, the normalization of the 
dipole cross section was fitted to data and two distinct 
initial conditions, inspired in the Golec Biernat-Wusthoff 
(GBW) [H and McLerran-Venugopalan (MV) [H mod- 
els, were considered. The predictions resulted to be al- 
most independent of the initial conditions and, besides, 
it was observed that it is impossible to describe the ex- 
perimental data using only the linear limit of the BK 
equation, which is equivalent to Balitsky-Fadin-Kuraev- 
Lipatov (BFKL) equation [s^. In next section we will 
compare the results of the RC BK approach to the ex- 



4 



perimental data on exclusive processes at DESY-HERA 
and present our predictions for the kinematical range of 
the future electron - proton collider [iB] ■ 



III. RESULTS AND DISCUSSIONS 

In what follows we calculate the exclusive observables 
using as input in our calculations the solution of the RC 
BK evolution equation. In particular, we make use of 
the public- use code available in [36j . In numerical calcu- 
lations we have considered the GBW initial condition for 
the evolution (we quote Ref. [Q] for details) and it was 
verified the MV initial condition gives cross section with 
overall normalization 10 — 15 % smaller and unchanged 
energy dependence. Furthermore, we compare the RC 
BK predictions with those from the non-forward satura- 
tion model of Ref. [l^ (hereafter MPS model), which 
captures the main features of the dependence on energy, 
virtual photon virtuality and momentum transfer t. In 
the MPS model, the elementary elastic amplitude for 
dipole interaction is given by. 



27ri?; e 



-Bit 



lAA(rQsat(a;,|t|),x), (9) 



with the asymptotic behaviors Qg^^.{x,A) ^ 
max((5§,A^) exp[— Aln(a:)]. Specifically, the t de- 
pendence of the saturation scale is parametrised as 



qL i^, \t\) 



?o(l + c|i|) 



A 



(10) 



in order to interpolate smoothly between the small and 
intermediate transfer regions. For the parameter B we 
use the value B = 3.754 GeV"^ [11]. Finally, the scal- 
ing function J\f is obtained from the forward saturation 
model [33|. 

Here, in order to take into account the skewedness cor- 
rection, in the limit that x' <C a; <C 1, the elastic differ- 
ential cross section should be multiplied by a factor Rg, 
given by , 38l | 



_ 22Ae+3 r(A, -f 5/2) 

r(Ae-t-4) ' 

din [A{x, Q^, A)] 



with A 



(11) 



d\n{l/x) 

which gives an important contribution mostly at large 
virtualities. In addition, we will take into account the 
correction for real part of the amplitude, using dispersion 
relations ReA/ ImA = tan(7rAe/2). In the MPS model, 
the skewedness correction is absorbed in the model pa- 
rameters and only real part of amplitude will be consid- 
ered. 

Let us start to compare the RC BK predictions to 
the available HERA data for exclusive vector meson {p, 
J/^" and T) photo and electroproduction. In Fig. [l]we 
present the predictions of the RC BK and MPS models 
for the diffractive p° vector meson production and com- 
pare it with the current experimental data from ZEUS 
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FIG. 3: (Color online) Energy dependence of the jp cross 
section for T photoproduction. Data from ZEUS and HI 
collaborations [iil. 14511. 



(left panel) and HI (right panel) Collaborations [39|, |40j . 
These measurements are interesting as they cover mo- 
menta scale that are in the transition region between per- 
turbative and nonperturbative physics, where saturation 
effects is expected to play an very important role. As 
the numerical RC BK solution there exists only for for- 
ward dipole-target amplitude we need an approximation 
to compute the non-forward amplitude. Here, we assume 
the usual exponential ansatz for the t-dependence which 
implies that the total cross-section is given by 



(12) 



1 


d(TT 




^ d(TL 






dt 


H 


dt 


t=0- 



Notice that values of the slope parameter By are not 
very accurately measured. We use the parametrisation 



Bv (Q' 



0.60 



14 



M2)0.26 



1 



(13) 



obtained from a fit to experimental data referred in Ref. 
[4l| . The uncertainty in this approximation can be larger 
than 20-30 % depending on the value. It is verified 
that the effective power Ae is similar for both RC BK 
(solid line curves) and MPS (long dashed curves) pre- 
dictions, with the deviation starting only at the higher 
values where the predictions differ by a factor 1.5. 
This can be a result of the similar small- a; behaviour for 
both models, where the effective power ranges from the 
soft Pomeron intercept Xe{Q^ = 0) ~ aip{Q) = 1.08 up 
to a hard QCD intercept Xe{Q^) — cNcag/Tr ~ 0.3 for 
large Q^. The data description is fairly good, with the 
main theoretical uncertainty associated to the choice of 
the light cone wavefunction (about a 15 % error). It was 
verified that the contribution of real part of amplitude 
and skewedness are very small for p production. 

In Fig. [5] we present the predictions of the RC BK 
model for the diffractive J/"^ production and compare 
with the ZEUS (left panel) and HI (right panel) data 
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FIG. 4: (Color online) Energy dependence of the DVCS cross 
section for different photon virtualities. Data from HI collab- 
oration [4611. 



[39l kOj . It is verified that the effective power Ae is simi- 
lar for both RC BK and MPS only in the photoproduc- 
tion case. The situation changes when the photon vir- 
tuality increases. The effective power for RC BK (solid 
line curves) is enhanced in in comparison with the 
non-forward saturation model (long dashed curves). The 
data description is reasonable since it is a parameter-free 
calculation and the uncertainties are similar as for p pro- 
duction. For J/^ production, the contribution of real 
part of amplitude increase by 10 % the overall normal- 
ization, while the skewedness have a 20 % effect. In the 
MPS model, as discussed before, the off-forward effects 
are absorbed in the parameters of model. The RC BK 
and MPS predictions differ by a factor 1.4 for large en- 
ergies. For sake of completeness, in Fig. [3] the results 
for T photoproduction is presented. The RC BK and 
MPS predictions are similar in the HERA energy range 
and differ by a factor 1.5 for large energies. It is known 
so far that the dipole approach underestimates the ex- 
perimental data for T and the reason is not completely 
clear [3, [l3]- However, the deviation concerns only to 
overall normalization, whereas the energy dependence is 
fairly described. The referred enhancement in the effec- 
tive power Ae is already evident in T photoproduction as 
the meson mass, my = 9.46 GeV, is a scale hard enough 
for deviations to be present. Skewedness is huge in the 
T case, giving a factor i?^ 1.3 in photoproduction. For 
this reason, we have included this effect in both mod- 



els. However, this is not enough to bring the theoretical 
results closer to experimental measurements. 

Finally, we analyse the DVCS cross section and com- 
pare it to the recent HI data f46l|. The cross sec- 
tions are presented as a function of W , for different val- 
ues of , in Fig. H) Here, the approximations con- 
cerning the final state particle are not present and the 
cross section suffers of less uncertainties. For the slope 
value, we take the experimental parametrization |4(| . 
B (Q2) ^ a[l-61og(QVQ2)]^ ^j^i^ ^ ^ 6.98±0.54 GeV^, 

6 = 0.12 ± 0.03 and Ql ^ 2 GeV^. The situation for 
DVCS is similar as for vector meson photoproduction, 
where the effective power Ae is similar for both RC BK 
and MPS for small virtualities and starts to change as 
grows. Skewedness is increasingly important for DVCS 
at high and it was introduced for RC BK model. For 
the MPS model this effect is absorbed in the its parame- 
ters as noticed before. The RC BK and MPS predictions 
are similar for the HERA energy range, describing the 
current data, and differ by a factor 1.2 for large energies. 

As a summary, we presented a systematic analysis of 
exclusive production in small-x deep inelastic scattering 
in terms of the non-linear QCD dynamics. This approach 
was performed using the recent calculation of the running 
coupling corrections to the BK equation. In this work 
we obtained the predictions for the exclusive production 
of vector mesons and DVCS and compare them to the 
available experimental results and the predictions of the 
MPS model. The main novelty of this work with respect 
to previous phenomenological analyses is the direct use 
of the running coupling BK equation to describe the en- 
ergy and virtuality dependences of exclusive processes at 
DESY-HERA. We find a fairly good agreement with ex- 
perimental data using a parameter-free calculation (pa- 
rameters are fixed from structure function F2 data). Our 
main result is that the RC BK evolution equation implies 
larger cross sections for exclusive processes than the phe- 
nomenological model proposed in (isj . Our predictions 
for both vector meson and DVCS production are rele- 
vant for the physics programs in the ongoing experiment 
LHeC and in the photoproduction processes in coherent 
proton - proton interactions at the LHC. 
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